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A new definition of coarse-grained quantities describing the dust flow in General Relativity is 
proposed. It assigns the coarse-grained expansion, shear and vorticity to finite-size comoving do- 
mains of fiuid in a covariant, coordinate-independent manner. The coarse-grained quantities are all 
quasi-local functionals, depending only on the geometry of the boundary of the considered domain. 
They can be thought of as relativistic generalizations of simple volume averages of local quantities 
in a fiat space. The procedure is based on the isometric embedding theorem for 5*^ surfaces and 
thus requires the boundary of the domain in question to have spherical topology and positive scalar 
curvature. We prove that in the limit of infinitesimally small volume the proposed quantities re- 
produce the local expansion, shear and vorticity. In case of irrotational flow we derive the time 
I evolution equation for the coarse-grained quantities and show that its structure is very similar to 

the evolution equation for their local counterparts. Additional terms appearing in it may serve as 
5— ( ' a measure of the backreacton of small-scale inhomogeneities of the flow on the large-scale motion 

^ ' of the fluid inside the domain and therefore the result may be interesting in the context of the cos- 

mological backreaction problem. We also consider the application of the proposed coarse-graining 
procedure to a number of known exact solutions of Einstein equations with dust and show that it 
(N ■ yields reasonable results. 
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I. INTRODUCTION 



' During the last few years there has been a sudden growth of interest in the problem of averaging of the Einstein 
equations of gravitation. The topic attracted the attention of the relativist community because the influence of 
small-scale inhomogeneities in the Universe upon its large-scale behavior, neglected by most cosmologists so far, may 
, possibly offer a simple and attractive explanation of the accelerated expansion of the Universe ( [H-Q , though see also 
' Q and 01 for a different view). The discussion seems inconclusive so far, but apart from the application to cosmology 
I the problem of averaging the Einstein equations is very interesting by itself, as it touches upon a fundamental question 
■ concerning the application of Einstein's equation to physical and astrophysical situations. Namely, when we apply 
^\ ' Einstein's theory to objects like galaxies or galaxy clusters, often treating the stars, black holes and other bodies of 
which they are composed as dust or fluid, we also effectively replace a complicated metric, describing the fine details 
^ ■ of spacetime in the vicinity of those objects, by a smooth one. When one realizes that, the fundamental questions of 
• ^ validity of such an approximation and corrections one should possibly include arise immediately j8|, lOj . 

Up to now the study of cosmological backreaction has mostly proceeded along two main lines of research. The first 
, is the study of several families of exact, non-homogeneous solutions of Einstein's equations filled with dust of fluid 
■ - - I (Lemaitre-Tolman-Bondi [lol - [l3| . Swiss cheese }14i - il6i] , Szekeres [12J etc.). The second one consists of attempts to 
derive a general expression for backreaction in a non-homogeneous metric, using perturbation expansion and/or a 
coarse-graining ( "averaging" ) scheme [l7l - l22| . In this paper we follow the second line and suggest a new approach to 
deriving an exact, non-perturbative formula for backreaction. 

The averaging scheme for the metric and matter fields is a crucial element of any backreaction study. It allows for 
smoothing the fine structure of a given spacetime and matter fields, pushing the entire configuration of fields towards a 
more symmetric and homogeneous one. In statistical physics the procedure of forgetting the fine, small-scale structure 
of a physical system is known as coarse-graining. J54| 

Clearly there are infinitely many ways to perform coarse-graining and the approach one chooses must be adapted 
to the physical system under consideration. In the context of cosmology, or broadly speaking General Relativity, 
several schemes have been proposed in the past. The most popular one is the Buchert's scheme [23j . in which one 
assigns scalar cosmological parameters - the volume expansion, the average scalar curvature and matter density - 
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to a comoving domain in 3+1 splitting, using its total volume and volume averaging. The resulting formalism is 
quite simple, but it suffers from a fundamental drawback - it only allows for coarse-graining of the scalar part of 
the evolution equations, completely ignoring the tensorial part. Zalaletdinov [1,[23|; proposed a framework based on 
a globally defined field of "bivectors" , connecting tangent spaces of distant points in spacetime. It is more general 
than Buchert's and allows for coarse-graining of the full Einstein equations, but it requires an earlier choice of the 
bivectors, which introduces certain unwanted arbitrariness into the results. It was later developed and applied to 
the cosmological context by Coley, Pelavas and Zalaletdinov [1^, H^, as well as Singh and Paranjape [Ij]. Other 
approaches have been proposed by Anastopoulos 28] , Carfora and Marzuoli 29] , based on the Ricci flow deformation, 
later developed in [30l - [32| . and Reiris [ssj . 

In this paper we will propose an alternative to Buchert's and other coarse-graining schemes. We will present a 
procedure which assigns the coarse-grained expansion, shear and vorticity to a finite-sized comoving domain of the 
fluid in a coordinate-independent manner. The proposed quantities are quasi-local, i.e. they are functionals of the 
geometry of the boundary of the coarse-graining domain. While quasi-local quantities are abundant in mathematical 
relativity [s^], they haven't been used too much in this particular context (though see [s^, [IB]). The shear and 
expansion are functionals of the induced metric on the boundary and its time derivative, while the definition of 
vorticity involves the shift vector field at the boundary. For computational simplicity we develop our formalism with 
pressureless dust as the only matter field present, but the construction carries over easily to more complicated fluids. 

The definitions proposed here are generalizations of the volume averages of shear, expansion and vorticity in 
non-relativistic, Newtonian theory of gravitating fluids. The basic idea is the following: we embed the boundary 
of a comoving domain isometrically into the three-dimensional Euclidean space and construct a fictitious, three- 
dimensional fluid velocity field on the image of the embedding. The velocity field is chosen in such a way that its 
'fictitious' flow induces the same infinitesimal metric deformation on the embedded surface as the 'true' dust flow does 
on the domain boundary in the spacetime. This fictitous velocity is then used to assign the coarse-grained expansion 
and shear via the standard, "Newtonian" volume average of the velocity gradient, transformed into a boundary 
integral. The vorticity on the other hand cannot be read out from the boundary metric only, so in order to define 
it we use a different approach, involving the pushforward of the ADM shift vector to the Euclidean space. In both 
cases the construction and the uniqueness of results rely on the existence and rigidity of the isometric embedding, 
guaranteed by the well-known theorem conjectured by Weyl and established by Lewy, Nirenberg, Heinz, Pogorelov, 
Alexandrov and Cohn-Vossen (see [s^, [ss'] and references therein) . The embedding theorem imposes two relatively 
mild restrictions on the coarse-graining domain, namely its boundary must be homeomorphic to a two-sphere and 
must have a positive curvature. The first one is quite harmless, as most interesting domains are ball-shaped anyway 
and thus have a topologically spherical boundary. The second should be satisfied if the domain is "round enough" , 
i.e. its shape is close to a geodesic ball or ellipsoid. 

The strength of the approach presented in this paper lies mainly in the fact that it is manifestly covariant, as 
each step can be defined in a geometric, coordinate-independent manner. Apart from the domain itself and the 3+1 
splitting of the spacetime, one does not have to make any additional, artificial choices which might influence the final 
result, unlike the approaches of Anastopoulos or Zalaletdinov. A clever use of the embedding theorem allows one to 
bypass the problems of covariant averaging of tensor fields and consequently the resulting decomposition of the fiow 
into the large-scale, coarse-grained part and local inhomogeneities is covariant. It is also worth mentioning that the 
expression for backreaction one finally obtains has the form of a surface integral divided by volume. This makes it 
relatively easy to estimate its magnitude given estimates on the magnitude of local inhomogeneities. 

The main purpose of coarse-graining over a domain is to find a best fitting homogeneous cosmological solution. 
In Buchert's original scheme the fitting is performed to a FLRW metric via three scalar parameters, one of them 
being the suitably defined expansion. The advantage of including the tensorial shear and vorticity is that it expands 
the class of homogeneous solutions one might want to fit to non-isotropic ones. The presented procedure is more 
mathematically involved than Buchert's, but still it seems to be a good compromise between mathematical elegance 
and computational simplicity, as it only requires solving two systems of well-known PDE's and evaluating a surface 
integral. 

The paper is organized as follows. In the next section we will briefly discuss the theory of Newtonian gravitating 
dust ("Newtonian cosmology"), the coarse-graining of its evolution equations and discuss the backreaction. The 
section is intended to be a short review and summary of the most important results, for an exhaustive treatment 
of the subject see [3^ 113] and j4l| - |43j . It is included to provide the motivation for the relativistic coarse-graining 
procedure. 

In the third section we will spell out the definitions of coarse-grained expansion, shear and vorticity for any spacetime 
filled with dust as well as the coarse-grained counterpart of the Fermi derivative operator for irrotational dust. In 
the fourth one we will prove that the coarse-grained quantities coincide with the local ones if the coarse-graining 
domain shrinks to a point, thus showing that the definitions are reasonable. In the fifth section we will derive the 
time evolution equation for the coarse-grained quantities for irrotational dust and show that it has a similar structure 
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to the local one. We will derive thereby an exact, non-perturbative expression for backreaction in a spacetime filled 
by non-rotating dust. The rotational case will be discussed in the next paper. 

Finally we will discuss briefly the application of the presented formalism to several known exact cosmological 
solution of Einstein's equations. 



II. NEWTONIAN COSMOLOGY 

Consider a cloud of dust in ordinary Euclidean space E^, interacting via Newtonian gravity. In Cartesian coordinates 
it is described by the local density function p{x°-,t), velocity field v°'{x^,t) and the Newtonian potential (f){x°',t). 
They satisfy a well-known system of PDE's 

dp „ dp dv°' 

A0 = 47rGp. (3) 

The first one is the equation of motion of the dust particles in the gravitational field, the second is the continuity 
equation and the third is the Poisson equation for the Newtonian potential. The left hand side of equations ([1} and 
([2]) is usually called the convective derivative and denoted by and respectively. 
Consider the solutions to (II|)-([3]) of the form 

v''{x'',t) = Q\{t)x'' 

p{x'',t) = p{t) 

Solutions of this type are homogeneous in the sense that the matter density is constant in space and the matter 

undergoes a homogeneous expansion, stretching and rotation. The ansatz plugged into the equations of motion it 
yields a system of non-linear ODE's 

Q\ = -Q\Q\-^\ (4) 

p - -pQ\ (5) 

<f\ = AnGp. (6) 

for one general matrix a symmetric matrix and a scalar function p. can be decomposed into the trace, 
called expansion, the symmetric tracesless part (shear) and the antisymmetric part (vorticity) 

Qab ~ -^SSab + (J{ab) + ^[ab] ■ 

Similar procedure applied to $a6 yields the scalar part $, directly related to p via ([5]), and the traceless part 
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b -I- ^(ab)- 



E(jf,(i) can be interpreted as the large-scale tidal field. It doesnt't have any evolution equations of its own, which 
makes the system (U])-® undcrdetermined. We can close it up if we assume the tidal field to vanish. 

The homogeneous solutions discussed here contain as a special case the Newtonian counterparts of FLRW solutions. 
It suffices to assume that the flow is rotation-free and homogeneous {Qab — H Sab) and equations (l4])-([6]) yield the 
standard Friedmann equations with dust matter. 



A. Coarse— graining of the Newtonian cosmology 

In order to calculate backreaction we need to split all relevant quantities into the large-scale, coarse-grained part 
and small-scale inhomogeneities. Fix a compact domain Gt G E'^, dragged by the dust particles. We assign coarse- 
grained quantities to Gt by simple volume averaging 



{Q\) ^ h I <6d'^ (7) 

IGt 



V 
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kab) = TP [ (b,abd^X (8) 
■/ Gt 



{P) - h I Pd'^c, (9) 



V 



Gt 



where V denotes the time-dependent volume of Gt- 

Note that by the virtue of the divergence theorem the first two averages are effectively surface integrals: 



{Q\) - ^ / v^nbda (10) 



X _ 1 
ab) - y 



(f>^a nb da. 

dGt 



They are therefore only sensitive to the values of w° and at the boundary of the coarse-graining domain. We can now 
split the velocity, density and the Newtonian potential into the large-scale part and the small-scale inhomogeneities: 

P ^ {p) + Sp- 

Now the problem of backreaction in Newtonian cosmology can be formulated as a precise mathematical question: 
how do the evolution equations for the coarse-grained quantities differ from those obtained with the assumption of 
complete homogeneity, i.e. (jj])-®? 

It turns out that despite the nonlinearity of the system (HJ-Q the former have the same structure as the latter 
ones 

{Q\)- = -{Q\){Q\)~{r,b) + B\ (11) 
{p) - -{p){Q\) (12) 

(^■%) - 47rG(p) (13) 
apart from just one additional term in the first equation: 

B\^V-^ [ iSv''6v''tnc-Sv''Sv'' cnb)da. (14) 

JdGt 

Clearly this is a backreaction term, as it describes the infiuence of the inhomogeneities upon the dynamics of the 
large-scale quantities. Equation was first derived in the special case of Q"^ = ^ ^"^b Ehlers and Buchert [s^. 

Expression (jl4[) for backreaction has a number of interesting features. First, note that it depends on velocity 
inhomogeneities, but not on the density inhomogeneities. Consequently the backreaction is exactly the same for a 
smooth matter distribution and a collection of point particles. The physical reason behind this surprising fact is the 
equivalence principle: in absence of pressure the motion of a particle is only affected by gravitational forces, which in 
turn does not depend on its mass. Therefore p drops out of equation ([T]) and Sp does not appear in the expression 
for backreaction. 

Secondly, Bab depends only on the values of 6v°' at dGt. Although (ITi|) appears to involve its derivatives in all 
directions, the normal derivative drops out because of antisymmetrization in b and c. 

Finally equation p4)) has the form of a surface integral divided by the enclosed volume. This is of course a pleasant 
consequence of the fact that the expressions for and {(j)^ab) are effectively surface integrals too. Thus (Q°f,)' 

depends only on the dynamics of particles at the boundary of the domain and any backreaction terms arising in that 
equation must be surface expressions too. 

This form makes equation (|14p a very convenient tool to provide bounds on backreaction on very large scales. 
Consider a ball-shaped domain of radius R. Integrating over the surface yields the factor of i?^, while by dividing by 
volume amounts to multiplying by R~^. If we can give absolute bounds on Sv"" and their derivatives on the spacetime, 
we can obtain even stronger bounds on the backreaction of the form of 

\B\\ < ^ 

for a positive constant C. 

Let us stress that the coarse-grained equations (jlip - (|13l) are exact, they do not involve any kind approximation 
or perturbation expansion. In fact, pT|) - (IT51) are just the rewriting of (IT])-® in terms of non-local variables ([T])-® 
and they hold universally for any solution. In particular, expression p4l) yields the correct value of backreaction no 
matter how much the solution considered differs from the corresponding homogeneous one. 
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III. COARSE-GRAINING IN GENERAL RELATIVITY 

The result for Newtonian backreaction, obtained in the previous section, suggests a new, promising way to approach 
the backreaction problem in the relativistic setting: one could derive the relativistic counterpart of equation (jlll) for 
the time derivative of a non-local, coarse-grained velocity gradient, involving additional terms coming from the metric 
and flow inhomogeneities. 

Consider a curved spacetime filled with gravitating dust, described by energy density p and four-velocity field . 
Define the gradient of 

Z^,v = Vixii^, (f5) 

which is completely orthogonal to 

Z^, u^" = Z,^u'' = 0. 

This makes effectively a three-dimensional tensor whose role is analogous to v"- from the previous section. It is 
conveniently decomposed into the local expansion, shear and vorticity 

which describe how an infinitesimal fluid or dust element is deformed and how it rotates during its motion. Z^^, 
satisfies an evolution equation along the worldlines of dust particles 

V„Z^, = -Z% Z", - R%,, u" u". (16) 

This equation is clearly analogous to (U) , since the contraction of the Riemann tensor on the right hand side can be 
decomposed into the trace, related to the Einstein tensor and thus to the local matter content, and the traceless Weyl 
tensor describing the tidal forces, just like 'f'"^ in (|4]). The trace of equation (IT6l) is the well-known Raychaudhuri 
equation. 

Let us try to push the analogy with Newtonian cosmology further and consider a finite fluid element traveling 
through spacetime. All particles enclosed within it make up a four-dimensional cylinder C in spacetime and particles 
contained in its boundary generate a three-dimensional tube dC (see fig. [IJ. The tube can be then foliated by 
previously chosen constant time slices. One could then propose a way of assigning the coarse-grained Z^i^ (and thus 
the expansion, shear and vorticity) to the constant time slices of C . The interior of these slices Ci would play the role 
of the comoving coarse-graining domain. We could then derive the time evolution equation for the coarse-grained 
{Z^J_y) and, if the proposed definition is appropriate, the evolution equation should have a similar form to with a 
bunch of additional terms. These new terms should arise from the inhomogeneities of the metric and the dust flow 
inside and on the boundary of the domain. They would describe how these small-scale inhomogeneities influence the 
large-scale motion of the coarse-graining domain as a whole. Like Bab in (IHp . they should vanish on a homogeneous 
FLRW solution but not in a general situation, which would justify referring to them as the backreaction terms. 

Obviously there are infinitely many possible ways to assign {Z^^-J) to a finite domain. We should however make sure 
that the prescription satisfies a couple of reasonability conditions. First, we the new definition should be consistent 
with p5|) . i.e. we should recover the local velocity gradient Vj/U^ if we shrink the tube under consideration to an 
infinitesimally small one. Second, a good coarse-graining prescription should be covariant in the following sense: 
apart from the choice of the tube itself and the 3+1 splitting of the spacetime, the result should not depend upon any 
externally introduced structure, including the coordinate system. This requirement is not just a question of aesthetics 
or mathematical pedantry, as any dependence of this kind introduces spurious degrees of freedom which must be 
separated out carefully to obtain meaningful, physical results. 

While scalar quantities like expansion are easy to coarse-grain in a covariant way by taking the volume average, no 
simple way to coarse-grain tensors exists in a general, curved spacetime. In fact, the covariance problem appears here 
at an even earlier stage: note that the coarse-grained vectors and tensors cannot live on their own, they need a vector 
space to be defined at all. A sensible prescription must begin by assigning a vector space to each coarse-graining 
domain, once again in a coordinate-independent way. 

We will now propose a relatively simple procedure of assigning {Z^y) to finite elements of the fluid, which can be 
proven to work under relatively mild conditions. We will prove that it satisfies the requirements stated above and 
derive the evolution equation for {Z^^y) in the special case of irrotational dust flow. More general case will be treated 
in the next papers. 

The main tool we are going to use is the isometric embedding theorem for S'^ surfaces: 
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FIG. 1: Four-dimensional cylinder C, generated by particles within a finite volume of the fluid, and its boundary dC, a 
three-dimensional tube of worldlines, both foliated by constant time slices. 

Theorem III.l (Isometric embedding theorem for S^) Given a compact, orientable surface S homeomorphic 
to 5^, with positive metric q whose scalar curvature R > 0. Then 

• there exists an isometric embedding 

into the 3~dimensional Euclidean space; 

• the embedding is unique up to rigid rotations, translations and reflexions. 

(see [131 fo'^ exhaustive treatment of the topic and the history of developments in this area; in particular the 
existence part is discussed in and the uniqueness in [i^). If we fix the orientation, the theorem states simply 
that a compact, two-dimensional surface of topology, satisfying the positive curvature condition, can be recognized 
as a surface in and moreover this can be done in only one way, up to moving the surface around or rotating it as 
a whole. Let us stress that the hypotheses of the isometric embedding theorem cannot be relaxed, the theorem is in 
general untrue if they are violated. Compact, orientable surfaces which satisfy them will be called admissible. 

The isometric embedding theorem allows one to recognize the boundary of the coarse-graining region dCt as a 
submanifold in E'^ as long as it remains admissible. dCt is at the same time the constant time section of the tube 
dC. We will show that one can read out (Z^^) just from the geometric data on dCt, completely ignoring the geometry 
inside the coarse-graining domain. While the idea of assigning a coarse-grained quantity to a finite volume in a 
quasi-local manner, i.e. using just the geometry of its two-dimensional boundary, may seem strange at first, note 
that this strategy worked remarkably well in the Newtonian case: ^ depends effectively only on the surface values 
of w", the behavior of the velocity field inside the domain simply does not influence the its value. 

A. Notation, 3+1 splitting, conventions regarding the indices 

In this paper we will mostly work in a comoving, though not necessary orthogonal coordinate system {t, j/'), i.e. in 
coordinates in which the following conditions 

u^y\^ = 

hold. A coordinate system of this kind should exist at least locally. Its constant time surfaces are spacelike and will 
be denoted by Et, while their orthonormal, future-pointing vector will be referred to as v"' . We can perform the 
standard ADM decomposition of the metric and curvature in this coordinate system, obtaining a positive definite 
3-metric hij, laps function N and shift vector iV*, orthogonal projection operator P"^ = 5°'^ + Ujj and the extrinsic 

curvature K^^ = ~P^a P^j ^ P ^he covariant derivative operator on Sj, adapted to hij, will be denoted by D. 
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FIG. 2: One-parameter sequence of embeddings and the fictitious Euclidean trajectory of a single particle 

In comoving coordinates the worldlines of all dust particles are the curves of constant spatial coordinates y'. We 
can therefore describe any tube dC of particle worldlines by equation of the type = where 9^ are two 

coordinates on the tube sections and is a set of three functions. The tube itself is now parameterized by the 

coordinate time t and 9^. 

The constant time sections of the tube will be denoted by dCt and the domain of St enclosed by them by Ct- The 
metric induced on dCt, expressed in {9^), is given by 

qAB{t,9'')^e^Ae,Bh^j- (17) 

As usual, the Greek indices (^, ^, ...) run from to 3 and will be used with geometric objects defined on the 
tangent space to the spacetime. The lower case Latin letters ... are assumed to run from 1 to 3 and will denote 
objects on the space tangent to constant time hypersurfaces S*. The lower case letters from the beginning of the 
alphabet (a,5, c, ...) also run from 1 to 3, but will be used with objects defined on the Euclidian space E'^. Finally 
the upper case Latin indices A, B, ... run form 1 to 2 and will be used with objects on two-dimensional boundaries of 
coarse-graining domains (constant time sections of the worldline tubes). 

B. The coarse— graining procedure 

Let dC denote a worldline tube whose constant time sections dCt are admissible in the sense of theorem IIII.ll 
Consider a family of time-dependent isometric embeddings 

ft : dCt >^ dDt C E^, 

where dDt is the image of ft. Let a;", a ~ 1...3 be the orthogonal, Cartesian coordinates in E'^. The embeddings can 
be described by equations 

x^^x''it,9^) (18) 
with three functions x"(t,0'^) satisfying the isometry condition 

qABit,9^) = SabX''^Ax'',B- (19) 

P9p is a non-linear partial differential equations system and while theorem llll. II guarantees the existence of solutions, 
they are in general impossible to find analytically. However numerical schemes for solving it have been developed in 
the context of black hole theory [i^ |43| . 

The isometric embeddings can be freely rotated and moved around in E"^ without violating 

X^it, 9^) ^ R\{t) x'it, 9^) + W^it). (20) 

For any choice of the isometric embeddings we can follow the image of each individual particle of dC , labeled by 
coordinates 0^, in the Euclidean space, obtaining thereby its fictitious trajectory in E'^. If the trajectory is regular 
enough, we can also assign to it a fictitious instantaneous velocity v°' at each time (fig. [2]). 

In fact, we do not need to create a whole family of embeddings to recover of each particle. For a given instance 
of time it is enough to have one embedding and the time derivative of the induced metric in comoving coordinates. 
This is guaranteed by the linearized version of the embedding theorem: 
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FIG. 3: Vector field, defined on a convex surface, inducing given metric change, is unique up to the restrictions of rotational 
and translational vector fields to that surface 

Theorem III. 2 Given a compact, orientable surface S homeomorphic to S*^, embedded isometrically into I? , whose 
scalar curvature R > 0, and a symmetric tensor field rAB on S . Then 

• there exists a vector field u", v'^{x'') G , defined on S C. , such that 

QAB = rAB 

when S dragged along v"" , 

• v"" is unique up to adding a vector field of the form 

Y" = Vf'^x^ + W'' (21) 
VLab = —^ba = const, W"" = const. 

In short, the theorem states that if we are given a single embedding of a surface and a time derivative of the metric 
in form of a symmetric tensor, we can always find the three-velocity field defined on the image of the surface, 
which induces the right metric deformation at the linear order when the surface is dragged along it. The vector field 
is unique up to a rotational and translational part, which constitute a six-dimensional vector space. Let us stress 
that v"^ is defined only at dDt, our procedure does not assign velocity vectors to points inside or outside the image of 
the embedding (fig. (H). 

Theorem IIII . 2 1 follows easily from lIIl!l1 if we consider a regular, one-parameter family of surfaces. The existence and 
the regularity of the family of embeddings imply the existence of the time derivative of the embeddings with respect 
to the time parameter. This derivative yields the vector in question. The uniqueness of the solution is guaranteed by 
the second part of theorem lIII.il (see also the "infinitesimal rigidity" theorem from fl^). 

The vector field v"" is related to qAB via the action of a differential operator V given by 

V[v^]^s = 2v\j,x\B)5ab (22) 
or equivalently, in a more geometric notation, 

V[v'']ab = 2v(A\B) + 2w„ Hab = qAB (23) 

where va denotes the lower-index projection of v to dDt, w„ is the outward normal part, the vertical line is the 
covariant derivative on dDt and Hab denotes its outward extrinsic curvature in E'^. Theorem IIII. 21 can in fact be 
understood as a statement about the operator V: it has a known, six-dimensional kernel of the form of (|21l) and a 
non-unique inverse . If we constraint the domain of V by six independent, linear conditions on w", the restriction 
of V will be one-to-one and invertible. 



FIG. 4: Canonical isometry g between the tangent space Ta;Ef and E'^, induced by the isometric embedding 



We now apply the introduced machinery to our problem. For a given t we obtain an embedding ft : dCt dDt C 
and a vector field u° defined on dDt. We can take the boundary expression ([TOl) to define the symmetric part of {Zab) 

{Z(ab)) = TT [ -"(a rib) da, (24) 

^0 J dDt 

where Vq is the Euclidean volume of the domain Dt enclosed by dDt- The key observation is that this symmetrized 
integral is insensitive to adding vector fields of type (|2ip . Thus, irrespective of remaining ambiguity of u", (j24p defines 
unambiguously the expansion and shear of the fiow, coarse-grained over the domain Ct, as tensors acting on the 
Euclidean space |55| . 

In order to coarse-grain the vorticity we follow a different strategy. First, note that in the comoving ADM coordi- 
nates the velocity gradient takes the form of 



(the minus sign in the last term is due to our convention concerning the order of indices in Zij ) . We have shown how 
to extract the coarse-grained version of the first two terms from the induced metric qab- The antisymmetric part, 
i.e. the vorticity, does not infiuence qab. Note however that it depends only on the shift vector N\ Therefore in 
order to define (Z^g^^), we will push iV' directly from the constant time slice St to E"^. 

The pushforward is done via a canonical isometry g between the tangent spaces T^^Et, x £ dCt, and E'^, induced 
by the embedding ft. It is defined as follows: we decompose any vector X G T^Tit into the normal and tangent part 
with respect to dCt. We then push the tangent part Xt G T^dCt by /t*, while for the normal part we demand that 
the product with the outward-pointing normal is the same in both spaces 

g{XT) = ft*XT 
g{XYna = X'n,. 

(see fig. S}. Note that g can be used to push any geometric object, tensor or vector, from T^Et to E^ and the other 
way round. It is straightforward to check that g is at each point x G dCt an isometry between T^Et and E'^. We can 
now define vorticity by 

{Z[ab]) = ^ I 9{N\anb}A<j. (25) 

1^0 J dDt 

It is possible and in fact convenient to use (pS)) to fix the rotational part of v"" by the condition 

/ V[an^da^ g{N)[anb]da. (26) 

JdDt JdDt 

With this choice the whole (Zab) is related to by equation 

(Zab) = ^ / VaUbda (27) 

^0 JdDt 
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and the analogy with ([7]) is complete. Moreover, if we choose the embedding ft at one time to, ([26)1 fixes them at all 
times up to irrelevant translations. These can be fixed as well, for example by demanding that the center of mass of 
dDf remains at the origin, but nothing in the presented the formalism depends on the postion of the surface within E'^ 
or on the constant part of Therefore the translational ambiguity of isometric embeddings is absolutely harmless. 
For the sake of brevity we introduce the following notation for surface integrals of type (l27l) : 

^0 JdDt 

We will also denote by [ ■ ] the unique inverse of V satisfying 

plus an irrelevant condition fixing the constant part W"". With this notation the definition of the coarse-grained Zij 
can be written down in a slightly more compact manner 



b] ' 



(28) 
(29) 



in which we have used the time derivative of ([T7)) to substitute ^Zi^^j-j C'^C"* b 'Iab- Note that the combination 



2 on Sf 



plays in our formalism the role of coarse-graining operator for Zi^j) , a symmetric tensors of rank 



Remark. Although the coarse-grained (Zab) is a generalization of the Newtonian volume average (Qab), it does not 
have a straightforward interpretation as a genuine average of the local over Ct . Equations ([^ and (^5)) can be 
written as volume integrals, just like (jlOp . if we extend the velocity field v"" to the interior of Dt in an arbitrary way, 
but since we did not define any mapping between the interiors of Ct and Dt, such an extension would not be related 
to Z^i, inside Ct in any obvious way. It seems therefore more appropriate to refer to {Zab) as a coarse-grained rather 
than average quantity. 



C. Preferred time derivative operator 

If we want to derive the equation for the time derivative of (Zab), we need yet another ingredient. The equation 
involves the time derivative of a tensor and there are obviously infinitely many different ways to differentiate it. In 
the infinitesimal counterpart of the discussed equation (ITS)) the role of the preferred time derivative is played by the 
covariant derivative along u'^, coinciding with the Fermi derivative along the particle trajectories. In this subsection 
we will introduce its generalization to the coarse-grained tensors in the case of irrotational fluid {uiij = 0), with an 
adapted 3+1 splitting for which it^ is orthogonal to Ej, i.e. A^* = 0, = 1. 

Let T""^- -^^ (t) be a time dependent but spatially constant field on the Euclidean space. In Cartesian coordinates 
in which condition holds we set the preferred derivative to be equal to the ordinary time derivative 

(9 

q^rpab... rjiab... 

cd..."^^ cd...- 

In a more general situation, in Cartesian coordinates x"^ in E"^ in which the image of embedding ft rotates one can 
prove that 

^T^'-cd... = §iT^\d...+W^T^-'---,,,,, + W\T-^---,,,,, + -.- + 



where 



Vo 



Wab^-TT / V[anh]dcr. (30) 



dDt 



measures the net rotation of dDt. 
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IV. SINGLE WORLDLINE LIMIT 



We will now justify the definitions of {Zah) and, if the dust is irrotational, the preferred time derivative given above 
by showing that in the limit of the tube shrinking to a single worldline we recover the standard 4-velocity gradient 
and Fermi derivative. 

Consider first a one parameter family of tubes given by 

y^^\i\e^). (31) 

where A is a positive parameter. In the limit of A — the tube shrinks to a single worldline = (see fig. (HI). At a 
given instance of time t = tg the induced metric qAB: expanded up to the subleading order terms in A, has the form 
of 

{t, e^) = \^ e.B (0) + A^ c.A e.B Kj,k (0) e (^^) + o(a^) (32) 

where hijifi) denotes the metric at point p with coordinates — 0. 

Let To denote the tangent space TpYito to the point p and let denote the canonical coordinates on it, related to 
by condition 

1 = W —— 

for any vector I ^ Tq. Tq is a vector space, naturally endowed with a flat, Euclidean metric inherited from Sf. In the 
linear but most likely non-orthogonal coordinates the metric has the form of hij{0). 

Since Tq is a three-dimensional Euclidean space, we may use it as the target space for embeddings ft^ for all A > 
(see once again fig. [5]). In the leading order of A the embeddings must have the form of 

where x"' can be expanded in terms of A as 

x'^it, e^,x) - xcio^) + y rv(o) e{o^) ^{0^) + o{x^), (33) 

as usual up to rotations and translations. In order to prove it it is sufficient to calculate the metric induced on the 
image of embedding ([55)1 . apply the formula for the Christoffel symbol and see that result agrees with l\'S'2\i up to 
O(A^); the second part of theorem IlII . 1 1 assures that ([55]) is unique up to rigid motions. 

Note that since Tq plays the role of the embedding space, we will refer to it by indices a, 6, . . . rather that i, j, . . ., 
which we will reserve tangent spaces at other points. 



A. Expansion and shear 



We will now calculate (Zj-^j-)), in the leading order in A, for the family of tubes defined above. Functions £_^{6^) 
do not depend on the coordinate time, but the metric hab{0) at p does. The time derivative of the induced metric is 
equal to 

qABiO"") = A2 habit, 0) r A e.B + 0{X^) (34) 
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in the leading order. The time derivative of the metric on the other hand is related to via 

kj ^-2N + 2Z?(, = 2Z(,,) . (35) 

in any comoving coordinate system. 

Consider now the vector field V in Tq given by 

V^w") = i /if,c(0) r^(0) + 0{\^) 

(O(A^) denotes here a term which is of the order of for fixed Euclidean coordinates w^, unlike in all other equations, 
where the limit A is taken for fixed boundary coordinates 6^). The restriction of to dDt, 

y-(^0A) ^ ^ /,^^(o) ^"{9^) + 0(A2), (36) 

induces the right change of metric in the leading order when the embedded surface is dragged along it: 

VIv'^Ub = QAB + 0(A3). 

This can be proven by substituting ((36|) to (|22]) and comparison with (p4|) . From the second part of theorem IIII.2[) 
we know that this is the only solution to up to Euclidean motions. We can plug this guessed vector field into 
(f24l) to obtain 

''O J dDt ^0 JDt 

= ^hab + OiX'') 

or 

(^(a;,)) = ^(a6)(0) + O(A2). 



B. Vorticity 

Calculating of the expansion and shear required only the retaining of the leading terms in A expansion. Vorticity 
requires including the next one as well, because, as we shall see, it is the subleading term in giN)"" which contributes 

to {Z[ab])- 

For convenience we introduce auxiliary coordinates 

^'' = 2/' + ^rV,(0)y^"/ 

in the neighborhood of p (z* — 0). They agree with ?/* in the first order at the origin 

(9z* 



y'=o 



but they are locally fiat there in the sense that hij^k{0) ~ 0. The coordinates on Tq are also compatible with , 
i.e. I = w°' for any vector / G Tq. 

The tube is described in the auxiliary coordinates by equation 



A2 
2 



This functional dependence is the same as in the equation (pS)) for isometric embedding in Tq: 

X"(0-^) = S'^(0^) + O(A3). 
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Consider the outward-pointing normal to dCt, given by 



nfc(r) = _ 

Vbi bm /i'™ 

{eijk stands for the totally antisymmetric symbol and the sign in the first equation depends on the orientation of 9"^). 
hk, expressed in z^, has the same functional dependence as the outward normal rio in E^, expressed in w'', 

up to the subleading order. This is a straightforward consequence of the fact that in the locally flat coordinates 
both the metric hij{z'') and the functions E'^{6'^) agree with hab{0) and x"(^^) up to the subleading order. Thus the 
pushforward of the shift vector agrees in these coordinates up to the subleading order as well: 

g{N)'' = N^iO) + 7V%(0) AC''(0'^) + 0{\^). 

Operator 7V[-] kills the constant part of S((iV)'* and decreases the order of other two terms by one, because the dtr is 
of the order of A^, while the volume Vq of the order of A^: 

= M [<7(iV)[a] ,] = %,6](0) + 0(A) 

Since the covariant derivative Da is torsion free, this means that 

= -D^aNt] + 0(A) = Z[ab]{0) + 0(A). 



C. Preferred time derivative operator 

Once again we assume in this subsection that the fluid does not rotate and A''* = 0. Let T^^-'-p^ (t) be a family 
of tensors at tangent spaces at points crossed by the geodesic = 0. We assume that the tensor field is orthogonal 
to u in each index. With this assumption T can be reconstructed unambiguously from its projection to subspaces 
TpSt c TpM. Such tensors are effectively three-dimensional objects. 

Consider the projection of V„T to TpE^ along u^^ 

where 11"^ — (5"^ +u°' up. In the ADM variables this derivative is equal to 

rpab... rvT) rpab... ryp 

~ ^ pd... ^ c~ ^ cp... ^ d~ 

We will now show that V coincides with the T) in the leading, zeroth order in A. 

Since T) has been defined in orthonormal coordinates, we first need to introduce them on Tq. Let 

be such coordinates at f = fo, i.e. let 

We can extend them to other times by solving the initial value problem 

L«,(to) = A"6 

for the ordinary differential equation 

= \hbc {L-'Y, S''" (37) 
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and taking = i°f,(t) uf' at all times. 

The antisymmetrized operator J\f vanishes when applied to the velocity field p6p . so in the introduced coordinates 
the matrix Wab defining V has the form of 

Wab^Q + 0{\). 

Any tensor field expressed in x"" can be converted back to w"" by the formula 

T'-' -ca... ^^"-^T^'-rs... {L-')\ {L-')\---L\L^,-.-. (38) 

If we apply ([551) to calculate PT"'' -^^ and make use of ([57)l . we obtain after simplification the equation for V in 
coordinates 

PT'"'-,,... = PT''^-,,... + 0(A). 



V. EVOLUTION EQUATION FOR {Zab) 



In this section we will derive the evolution equation for {Zab) and compare it with the evolution equation for Zij. 
We will only be concerned with irrotational dust, leaving the case of non-vanishing rotation to the next paper. 

In a spacetime filled with irrotational dust we can introduce comoving and orthogonal coordinates, in which N = 1 
and A^* = 0. Let Ct be a fixed, comoving domain, whose boundary is admissible at least for some time. We can 
introduce a one-parameter family of isometric embeddings satisfying and assign the coarse-grained shear and 
expansion to it, encoded in {Zab) = {Z(ab))- This allows for decomposition of the local Zab, pushed forward to 
dDt C E'^ via g, and the velocity field w° into the the coarse-grained, large-scale part and local inhomogeneities 



Zab 



{Z\)x^ 



- 5 Zab 



(39) 
(40) 



both defined on dDt only. 

We now pass to evaluating 'D{Zab)- We first note two identities concerning the time derivatives of operators M[-\ 
and V associated with dDt- The first one, 



^N[Xa], = -{Z^,}M[Xa],+M 



dt 



DXa 



dt 

+ M[v^-^,Xa].~M[v\,Xa]^ 



(41) 



is valid for any vector field X"^ defined on dDt for some time interval. Differential symbol ^ stands here for the 



substantial derivative ^ 



,,a d 



in E along the velocity v , it should not be confused with the covariant derivative 



on Ej. The second identity concerns the time derivative oiV ^ [ • ]. Let —V^ [fABT, then 



DY, 



dt 



- = V-'[fAB]a^V-'[2v\^V,^B\^ 



^b\cY[a] 



b] 



Wa 



(42) 



{Wa denotes, once again, an irrelevant constant vector, whose value depends on the condition we imposed on [ • ]). 
Both can be proven easily by taking the definitions of N and V^^ , differentiating them with respect to time and 
rearranging the terms. 

By applying (|4T|) to the velocity field itself, we get 



^M[Va], = -U[v\N[vb] 



M\v\v. 



DVg 

dt 



,b"<^ib- 

The convective derivative of v"- can be related to the second derivative of the metric qab using (|42 

'DVaJ 

b 



(43) 



dt 



b] 



(44) 
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QAB is in turn related to the time derivative of Zij 

QAB = 2Z^j y\B 
in ADM variables {t,y^). From ((T6| and (|35|) we evaluate this derivative as 



Zik Z 



R 



iOjO 



(45) 



and thus, putting (gS]), dH]), ^ and (gll) together, 
5 



J\fp- 



ah 



'^RiOjO C,A C"',_B 



a6 



The left hand side of this equation is just the preferred time derivative V{Zab): because Na = and TV [v[a] ^ = 0. It 
now suffices to plug in the decomposition (|39 )l -(|40 |) into the right hand side and rearrange the terms to obtain finally 
the evolution equation 



T^{Zab) — —{Zac) {Z'^b) ^ {RaObo) + Bab + Bab, 



(46) 



where 



{Rai 



ObO) 



'^RiOjO C,A 



{ab) 



is obviously the coarse-graining of the contraction of the Riemann tensor with u^, considered as a symmetric tensor 
on Et (compare with ([SS)) '). 

This is clearly the coarse-grained counterpart of equation (jl6p . It includes two additional backreaction terms; the 
first one is the symmetrized version of (jl4p . and we will call it the Newtonian backreaction 



Bab = ^f[Sv'^Jv(^a]^,•^ - A/" 



The second one is entirely relativistic and has a more complicated structure 



Ba 



Mr- 



ah 



[2 {5Z,, 5Z^, x\a x\b - Sv',A Sv',b Sa)] 



Just like in the Newtonian case, both backreaction terms are surface integrals divided by volume. In contrast to (|14p . 
however, Bab involves linear terms in perturbations. 



VI. APPLICATION TO EXACT SOLUTIONS 



We will now discuss briefly the application of the coarse-graining procedure to well-known exact solution of the 
Einstein equations with dust. In particular, we will show that in case of metrics where matter undergoes homogeneous 
expansion or shear the coarse-grained quantities behave in a predictable way, i.e. they are equal to their spatially 
constant, local counterparts irrespective of the domain under consideration. 

Consider first a FLRW metric with dust, closed, open or flat, in standard 3+1 decomposition 

g = -dt^ + a{tf h{t) (47) 
h = ^"f ^ + (d6'2 + sin^ 6 Aip'^) . 

We will show that for any domain Ct C St, admissible in the sense of theorem IIII.ll the coarse-grained velocity 
gradient consists only of the scalar part, equal to 3 times the Hubble parameter 



{Zab) ^ - Sab ^ H{t) Sab- 

a 



(48) 
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The vanishing of the vorticity fohows obviously from the vanishing of shift in (j47p . The metric induced on the 
boundary of any domain in undergoes a homogeneous expansion with time 

qAB(tj = qAB{to) 



whose momentary expansion rate is related to the Hubble parameter by 

qAB{t) = —qAB{t) = 2H{t)qAB- 
a 

It is straightforward to see that the same change of metric is induced on the image of the surface dDt in E'^ by a 
homogeneous expansion along the vector field 

V = H{t) x", 

restricted to dDt . (|48p follows now from ([M]) if we use the divergence theorem to turn it into a volume integral. In a 
similar fashion one can prove that in Kasner solutions [49| 

{Zab) = CFab + - 5ab, 

where 9 and aab are equal to the local expansion and shear. 

It is well-known that it is possible to replace one or more non-intersecting balls within a dust FLRW solution by 
balls excised from an appropriately chosen LTB, Schwarzschild or Szekeres solution (so called Swiss cheese models, 
see [13, mi, [l^l)- In this case {Zab) also consists only of expansion, related to the Hubble constant of the exterior 
FLRW solution as in p8)) . as long as the boundary dCt does not pass through any of the excised halls (though they 
can be present in any number in the interior of the domain). This follows from the result above and from the fact 
that the presented coarse-grained quantities are quasi-local and thus insensitive to the changes of the metric inside 
or outside dCt- Thus in swiss cheese models our procedure recovers for appropriately chosen but generic domains the 
cosmological Hubble parameter of the background homogeneous solutions. 

The only metric in which dust undergoes homogeneous, rigid rotation, without expansion or shear, is the well- 
known Godel solution 53] with cosmological constant. As one might expect, in this spacetime only the antisymmetric 
part of {Zab) does not vanish, though its exact value depends on the 3+1 splitting. 

A more detailed discussion of the application of the coarse-graining procedure to these and other exact solutions 
we be given in subsequent papers. 



VII. CONCLUDING REMARKS AND ACKNOWLEDGEMENTS 



In the subsequent papers we will derive the expression for backreaction for rotating dust and discuss more thoroughly 
its application to exact cosmological solutions. 
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